16867  BOSQUE  DRIVE 
ENCINO,  CAUF.  91436 

(2131  981-7683 


!i£{rasT/8o- 

Report  1-77 


An  Advanced  Automatic  Test  Generation 
System  for  Digital  Circuits^  A — 


Submitted  to 


Department  of  the  Navy 
Office  of  Naval  Research 
Arlington,  Virginia  22217 


Project  Director 


Melvin  A./Brouer 


NR  <*8-631 


This  work  was  carried  out  under 


PISTRIBimON_STATFWrMT 
Approved  for 


public  release; 
Distribution  U .limited 


TABLE  OF  CONTENTS 


ABSTRACT 


INTRODUCTION 


PREPROCESSING 


2.1  Rate  Analysis  . . 

2.2  Cost  Assignments 


TIMING 


3.1  Basie  Definitions  on  Time 43 

3.2  Problem/Solut ion  Space 30 

3*3  Dynamic  Assignment  of  Biases  and  Periods 34 

3.4  Data  Structures  A ......  38 

TEST  GENERATION  ALGORITHM  69 

4.1  Data  Structure 69 

4.2  Functional  Element  Routines  70 

4.3  Implication  Routine  74 

4.4  D- Drive  Routine  73 

4.3  Justification  Routine  ......  76 

4.6  Initialization  Routine 77 

4.7  Consistency  Routine 77 

4.8  Definition  of  D- Frontier 77 

4.9  Initialization  Module  77 

4.10  Preprocessing  79 

4.11  Unsolvable  States  •••••  ...  80 

4.12  Solvable  States  8l 

4.13  Backtrack  State  Event 82 

4.14  Flow  Charts  for  STIM  GEN 83 

FAULT  SIMULATION 92 

5.1  Introduction  92 

5.2  Simulation  System  Structure 93 

3*3  Data  Structure  for  Circuit 96 


P| 


r 


Bection 


Simulator  Data  Structure  97 

5«5  Simulator  Operation  ..••••••  . . . . . 101 

5.6  Fault  Insertion  ......  105 

$•7  Initialisation ( 106 

5*6  Data  for  Primary  Inputs 106 

5.9  Input  Skew  Buffering  .......  106 

5*10  Miscellaneous  Iteas  .....  106 


w 


1 


ABSTRACT 


This  1s  the  third  and  final  report  documenting  oar  research  into 
the  design  of  a new  and  powerful  automatic  teat  generation  system  for 
digital  logic,  called  TEST/80.  This  system  is  designed  around  the  follow- 
ing six  concepts. 

1.  A powerful  circuit  preprocessing  analysis  which  leads  to  greater 
efficiency  during  stimulus  generation. 

2.  An  effective  initialization  algorithm. 

3«  The  use  of  time  frames,  phases  and  periods  so  that  asynchron- 
ous circuits  can  be  accurately  processed  during  stimulus 
generation. 

V.  The  use  of  functional  level  primitives  so  that  complex  circuits 
including  shift  registers,  counters  and  RAM's,  can  now  be 
effectively  processed. 

5»  The  use  of  a stimulus  generation  algorithm  which  incorporates 
the  concepts  in  1 - 4. 

6.  The  use  of  a functional  level  concurrent  fault  simulator,  used 

4 / y 

to  grade  a test  and  produce  a fault  dictionary. 

This  system  is  Intended  to  replace  existing  test  generation  systems^ 
used  by  the  United  States  Navy,  and  meet  the  Navy's  needs  for  automatic 
test  generation  during  the  next  six  to  ten  years.  At  that  time  it  is 
believed  that  circuit  densities  will  have  made  the  use  of  such  a program 
infeasible.  The  report  consists  of  the  results  of  our  research  into  the 
solution  of  the  six  concepts  previously  presented,  and  does  not  represent 
a detailed  design  of  TEST/80. 
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1.  UTROEUCTION 


This  report  documents  our  results  on  the  conceptual  design  of  a 
very  powerful  and  efficient  automatic  test  generation  program  (ATGP) 
software  system  for  digital  circuits.  The  general  structure  of  such  a 
system  is  shown  in  Figure  1.1.  This  structure  is  quite  universal,  and  is 
used  by  most  existing  systems  such  as  LOGOS,  LASAR,  FAIRTEST,  L- LOGIC, 

^ and  INDICATES.  The  input  to  these  systems  is  a manually  generated  circuit 

model.  The  model  is  then  preprocessed,  where  syntax  errors  are  detected, 
library  models  are  substituted  for  calls,  and  fault  collapsing  is  carried 
out,  etc. 

In  the  STIM  GEN  module  test  vectors  are  generated;  this  can  be  done 
either  manually,  by  random  technique,  or  via  some  heuristic  or  algorithmic 
I procedure . 

In  the  simulation  module  the  test  just  generated  is  evaluated,  i.e., 
the  set  of  an  faults  detected  is  determined.  If  the  percent  detection 
1 exceeds  some  bound,  then  a dictionary  can  be  generated  and  a test  tape 

produced.  If  the  percent  fault  detection  is  too  low,  more  tests  are 
generated. 

There  are  numerous  reasons  why  existing  systems  fail  to  be  effect- 
ive for  present  day  circuits.  These  problems  will  be  even  more  severe  as 
more  and  more  LB  I devices  continue  to  be  used. 

Same  of  these  problems  are: 

1.  Existing  ATGP  systems  fail  to  initialize  circuits  which  are 

*1  lnltiallzable. 

2 . Existing  ATGP  systems  cannot  effectively  process  circuits  con- 
taining RAM's  and  large  counters  and  shift  registers.  Th^.s 
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problem  manifests  itself  by:  a)  not  generating  tests  for 

faults , and  simultaneously  b)  using  up  a tremendous  amount  of 
CPU  time  in  searching  for  a test. 

3.  Existing  ATGP  systems  cannot  effectively  handle  asynchronous 
circuits. 

U.  Some  existing  ATGP  systems  use  low-level  logic  primitives  and 
this,  together  with  3 above,  require  complex  modeling  and  work- 
arounds in  order  to  generate  the  input  circuit  model. 

TEST/80  has  been  designed  to  eliminate  and/or  drastically  reduce  the 
problems  listed  above.  The  basic  concepts  employed  in  TEST/80  are: 

1.  Functional  level  modeling. 

2.  The  use  of  time  frames,  phases  and  periods  to  handle  the  genera- 
tion of  tests  for  asynchronous  circuits. 

3.  The  use  of  a powerful  initialization  process  so  that  complex 
circuits  can  be  automatically  initialized. 

H,  The  use  of  new  sophisticated  preprocessing  concepts  to  reduce 
CPU  time  during  stimulus  generation. 

5.  The  use  of  a functional-level  concurrent  - driven  fault 
simulator. 

Preliminary  hand  calculations  indicate  that  for  complex  circuits  the 
potential  performance  of  this  system  will  be  significantly  greater  than 
existing  systems.  The  Improvement  in  performance  is  a function  of  the 
structure  of  the  circuit  to  be  tested,  and  hence,  no  absolute  number  can 
be  stated.  For  combinational  logic,  our  system  will  actually  be  a little 
slower  than  existing  systems  because  of  the  overhead  in  our  preprocessing 

. 
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and  complex  data  structures.  For  circuits  having  counters  and  complex 
decoding  circuitry,  our  system  can  reduce  a 5-hour  (1108)  test  generation 
problem  to  a 30-minute  (ll08)  test  generation  problem. 

This  report  is  the  third  and  final  report  on  TEST/80  to  be  published 
under  the  present  contract.  The  previous  two  reports  are: 

1.  TEST/80,  Report  No.  1-76:  "Functional  Level  Modeling  of  Complex 

Elements  in  TEST/80." 

2.  TEST/80,  Report  No.  2-76:  "initialization  of  Digital  Logic 

Circuits." 

The  first  report  presents  our  results  for  developing  functional  level 
models  for  complex  sequential  devices  to  be  used  by  our  stimulus  generation 
and  simulation  modules . 

The  use  of  these  high-level  primitive  models  reduces  CPU  time  during 
test  generation.  If  one  is  interested  in  an  internal  IC  fault,  the  module 
containing  the  fault  is  modeled  at  a lower  level  such  that  the  fault  is 
external  to  the  primitive  elements  in  the  model.  The  remaining  elements 
in  the  circuit  can  still  be  modeled  at  the  higher  level. 

The  second  report  deals  with  the  development  of  new  and  powerful 
algorithms  for  carrying  out  automatic  initialization  of  circuits.  Unfor- 
tunately, for  some  complex  circuits  these  procedures  may  not  be  computa- 
tionally feasible.  We,  in  these  instances,  plan  to  use  classical  initial- 
ization techniques  based  upon  three  (0,l,u)- valued  logic. 

The  present  report  deals  with  the  following  four  major  areas: 


( 


**  Preprocessing 

Rate  analysis 
Cost  analysis 
Timing 

**  Test  Generation  Algorithm 

**  Functional  level  concurrent  fault  simulation. 


We  only  deal  with  the  "research  aspects"  of  the  TEST/80  system,  i.e., 
with  those  problems  investigated  and  solved  during  the  course  of  our 
research. 
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We  have  not  included  any  discussion  of  the  classical  components  of 
an  ATGP  system  whose  solutions  are  well  known,  such  as  fault  collapsing, 

dictionary  construction,  etc. 
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2.  PREPROCESSING 
2 .1  Rate  Analysis 

Rate  analysis  is  a preprocessing  technique  in  which  certain  lines  in 
a circuit  are  assigned  labels,  called  rates.  These  labels  indicate  the 
maximum  rate  at  which  a line  can  change  values.  These  rates  are  used 
during  test  generation  in  one  of  two  ways. 

1)  Logic  values  can  be  assigned  to  lines  consistent  with  the 
maximum  rate,  and 

2 ) sequences  which  cannot  be  justified  can  be  identified. 

» 

Let  a and  b be  rates . Then  the  following  are  rates : 

1)  0 and  1 

2)  ab  (concatenation) 

3)  an  ■ a.a  •••  a 

n times 

M (a) 

5)  a*  ■ a *e  • • • (indefinite) 

6)  r (a  random  bit) 

7)  rn  (a  randan  sequence  of  0's  and  l's). 

Examples : 

a) 

b) 

c) 

-6- 


(01)*  - 0 1010  1 »• 

(oW  - 001110011100111  ••• 

(o2 (oi2 )2)*  - 0001101100011011  ••• 


0 


Normally,  when  no  confusion  will  occur,  we  do  not  indicate  the  *. 

If  rate  a has  length  / then  we  denote  this  fact  by  a[i]. 

Examples : 

a)  0314[7] 

b)  (02(012)3[9])2  [22]. 

A rate  of  the  form  e e is  of  order  1,  where  e e {0,1} . If 
rates  a and  b sure  of  order  i and  j respectively,  then  rate  ab 
is  of  order  i + j.  A rate  indicates  the  maximum  rate  at  which  a signal 
on  a line  can  change  values.  We  desire  to  approximate  these  maximum  rates 
by  simple  rates  of  low  order,  say,  of  order  1,  2,  or  3. 

If  A is  a signal  line,  then  the  rate  associated  with  the  line  is 

denoted  by  A . 

r 

Expansion : 

Let  b be  a rate.  Then  the  expansion  of  b,  denoted  by  E(b),  is 
also  a rate,  and  is  obtained  by  replacing  any  element  x € [0,1]  in  b by 
xn  for  an  arbitrary  n > 0.  Note  that  the  expansion  operator  can  be  used 

repeatedly. 

Example:  E(o2l)  - 0514  . 

We  say  that  b is  less  than  E(b),  denoted  by  b < E(b),  for  all 
expansions  of  b. 
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Compatible : 

Let  a and  b be  two  rates.  Then  a is  compatible  with  b if 
there  exists  an  expansion  of  b,  denoted  by  E(b),  such  that  a - E(b). 

If  a is  compatible  with  b,  we  write  b < a. 

Examples : 

a)  » »0  0 1 0 1 is  not  compatible  with  b ■ 021^. 

2 3 

b)  a * 0 1 is  compatible  with  b * 0 0 1 0 1.  This  follows  be- 
cause b can  be  expanded  as  follows: 

b - 

Given  a rate  a,  then  b is  the  greatest  lower  bound  of  order  i 
of  a if  ' 

a)  b is  of  degree  i 

b)  a is  compatible  with  b 

c)  there  exists  no  other  rate  c of  degree  i such  that  b is 
compatible  with  c,  i.e.,  there  exists  no  c such  that 

c < b < a • 

Consider  the  rate  A shown  below  for  line  A. 

r 

A - (OCWllOOOlllOOOOlllOOOllOOOll)  . 

Since  the  minimum  number  of  consecutive  l's  and  0's  is  2 and  3, 


00101 
I N^‘ 


0011100111  - (0  1*7 


,2,3x2 
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respectively,  the  greatest  lower  bound  for  Ay  is 


I If  one  tries  to  Justify  the  sequence  A ■ 0 0 1 0 1,  since 

E(A^)  / A for  any  expansion,  A cannot  be  justified. 

i 

The  maximum  rates  are  01  and  10,  and  the  minimum  rates  are  0 

* 

and  1 , namely  a constant.  We  assume  all  primary  imputs  have  maximum 
rate,  unless  specified  otherwise. 

The  complement  of  a rate  a is  obtained  by  changing  all  0's  to 
l’s  and  l's  to  0's,  and  is  denoted  by  a. 

The  AND  and  OR  of  two  rates  a and  b are  denoted  by  a*b 
and  a + b,  respectively,  and  are  also  rates. 


Rules  for  Propagation  of  Rates  Through  Gates 
1.  AND 


If  A ■ 1 then  C ■ B . 

r t t 

If  A * 0 then  C • 0 . 

T T 


i 
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If  A ■ 0 then  C - B 
r r r 


If  A - 1 then  C - 1 

T T 


3.  Inverter 


B 


> r-<\) 


Procedure  AND:  Computing  - Xr  • Yy 


Let 


and 


where 


T “l  -"2  n3 
X.  «x  x x ••• 

r 

“l  -®2 

Yr  - y y y5  ••• 

i1!  - “l  * 


Find  a prefix  of  Yr  of  length  n^.*  If  x - 0 then  output  0 , else 

°1 

output  the  prefix.  Delete  x from  Xr  and  the  prefix  from  Yy.  If 
processing  a term  of  the  form  ( •••  )n,  then  take  the  result  produced, 
say  p,  and  replace  it  by  (p)n. 

Return  to  the  new  Xy  and  Yp  rates,  relabel  if  necessary,  and 
repeat  the  procedure  until  both  have  been  processed  to  the  end.  In  order 
to  generate  low  order  rates,  one  can  restrict  this  procedure  to  the  case 


where 


n^  - or  m^+n^  or  m^  + st. 


+ n. 


, • • • 
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Example : Given 


_ ,,16,16  . v .2  6 

X - 0 1 and  Y ■ 0 1 

r r 


compute 


Z - 


X . Y 
r r 


1.  Ol6[l6] 


((02l6)[8])2[l6] 


.16 


2.  X^  - I16  - 116[16] 
r 

Y *=  0216  - (0216[8])2[16] 
r 


(o2i 6f 


zr  - 016(0216)2  . 


Example : 


oV6 

2 6 
0 1 


0l6(02l6)2  

01  — y 


l6(02(0l)3)2 


0l6(02(0r)5)2 


The  output  can  be  approximated  by  the  first-order  rate  O1^  r^. 


Creation  of  Rates 

Rates  are  usually  created  at  the  output  of  sequential  devices,  such 
as  flip-flops,  counters,  shift  registers,  etc*  Normally,  ve  select  the 
input  to  the  device  to  have  rates  such  that  the  outputs  are  maximum. 


Consider  a device  triggered  by  a positive  edge  on  the  clock  line. 

“l  °2 

We  denote  this  condition  by  Cp(t ).  If  (C^)r  ■ 0 1 , then  the 

rate  of  the  positive  edge  is 


«C  (t))r  ■ 0n  ln  where  n ■ + ng 


is  the  period  of  the  clock  pulse. 


JX  edge  triggered  flip-flop: 


1.  Maximum  rate  configuration 


1 — J 

“l  “2 
0 x 1 - C 


1 Hk 


onin 


in  on 


•“l+  "2 


Note:  For  n^  ■ n^  ■ 1,  n ■ 2 and  ■ 0 1 . This  illustrates 
the  fact  that  this  configuration  represents  a divide  by  2 counter. 


2.  General  configuration: 


"l  °2 
1 * 

oV* 

0-5  l"S 


..  oa  1^ 
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1)  a,b  > n j i.e.,  the  output  cannot  change  faster  than 

the  input  period 

2)  a > n^,  i.e.,  as  long  as  J - 0 the  F/F  cannot  be 

set 

5)  b > n^,  i.e.,  as  long  as  K » 0 the  F/F  cannot  be 

reset. 


3«  Special  case  (JK  F/F  acting  as  a delay  flip/flop) : 


r 


1.  Count  mode 


E - 1 
r 

(C  ) 

a 

oni° 

p r 

■ 

”l  + °2 

Vr 

■ 

0*1* 

m 

o2*  l2* 

m 

0^1* 

m 

08q  x8q 

c 

m 

0l6q-m  xr 

2.  Maximum  count  mode 

Set  n • m » 1,  q - 2. 


5.  Special  case 


“l  "2 
oi 

"l  + V "2  -2  * 


Then 

<Vr  - O”1" 

<s>,  - °2“ i2n 

etc. 
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k.  Parallel  load  mode : (i  f {A,B,C,D)) 

(Q. ) ■ L • Z vhere  Z la  the  rate  of  Q.  when  operated 

x r r r r i 

as  a delay  flip-flop  with  input  and  clock  C^. 

Shift  register: 


(Qi  )r  ■ Er  * zr  where  Zy  is  the  rate  of  a delay  flip-flop 
having  input  D and  clock  C^. 

Rate  Assignment  Procedure 

The  following  procedure  assigns  rates  to  the  lines  in  a circuit: 

t 

1.  Sequentially  process  clocked  elements. 

a 

2.  Given  an  element,  for  those  inputs  not  yet  assigned  a rate, 
assign  rates  to  achieve  maximum  output  rates. 

3.  Compute  output  rate  of  device. 

4.  Propagate  rates  through  combinational  logic  as  long  as  lower  order 
rates  can  be  preserved. 

5.  Return  to  the  beginning  and  process  the  rest  of  the  devices. 

If  the  input  to  a sequential  device  is  given  a new  rate,  reprocess 
the  element.  * 
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lcatlons  of  Rate  Analysis 

The  major  application  of  rate  analysis  is  the  identification  of 
allowable  sequences  on  lines. 

Example: 


Assume  E,  Cp  and  J are  driven  by  large  combinational  circuits 
and  hence  their  rates  are  not  known.  Let  K be  a primary  input. 

We  assign  rates  as  follows:  For  the  counter,  the  maximum  rate 

inputs  are  * 1 and  (Cp)r  - 01 . This  gives  us  Cy  = 031  1.  For 

the  JK  flip-flop,  the  maximum  rate  inputs  are  J ■ K ■ 1,  and  since 

(C  ) - 031  1,  we  have  that  Q - 032  l32 . 
p r r 

Assume  that  during  the  generation  we  desire  to  test  for  the  fault 
A s-a-0.  To  solve  this  problem  we  have  the  following  subproblems: 

1.  Set  A - 1 

2.  To  solve  1 we  have  subproblems: 

a)  t:  C - t (C  (t)  - 1,C  (t  - l)  - 0) 

P P P 

b)  t-1:  J - K - 1,  Q - 0,  Cp  - 0 . 

3.  To  solve  the  problem  Q ■ 0 we  have: 


© 


c)  t*  < t-1:  R - 0 


It.  By  Implication,  we  see  that  R ■ 0 also  resets  the  counter  and 
we  have  C ■ 0. 

In  summary  we  have, 

Time:  t'  •••  (t-l)  t 

Line  C:  0 •••  0 1 

But  rate  analysis  implies  that  « 0^  1.  Therefore,  for  t'  ■ 1,  we 
have  (t-l)  • 31  and  t ■ 32.  Therefore  a 32  clock  time  test  is 
required. 

Example: 

This  example  shows  how  our  functional  modeling  language  interacts 
with  rate  analysis  to  give  information  for  test  construction. 

1 
0 

(T  ■ translator ) 

To  set  A ■ 1 we  require: 

RfH14  U)4  (reset  and  count  up  4 times)  (l) 

' U (E  - 1 , C - t ) 
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r 

H ^ » (E  - 0 (impossible);  E ■ 1,  E ■ 01  } 

p 

R {D  - 0} 

Substituting  back  into  (l),  we  get  for  the  line  C : 

P 

C :R(h“  U)U  - x((0l)nt)4 

P 

Since  t * (01  ) and  by  rate  analysis  we  have  that  (C  = 0^  1,  we 
obtain  (0  l)n  = 0^°  and  hence,  Cp  = x 0^°  0 1. 

Another  application  of  rate  analysis  is  the  following:  Assume  that 

in  order  to  get  a test  for  a fault  we  assign  line  A the  sequence  A. 

Now,  if  A is  not  compatible  with  Ay,  then  A cannot  be  Justified  and 
t the  test  generation  procedure  should  backtrack  and  take  another  choice. 

2.2  Cost  Assignments 

In  preprocessing  a circuit  we  assign  to  each  line  A in  the  circuit 
three  costs,  namely, 

, cA  - the  cost  of  setting  line  A to  a 1 

cA  - the  cost  of  setting  line  A to  a 0 
and 

| dA  - the  cost  of  driving  a D (5)  on  line  A to  a primary 

output. 

These  costs  are  used  by  our  test  generation  algorithm  during  D- drive 

I 

and  line  Justification.  In  this  section  we  discuss  how  the  preprocessor 
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assigns  these  costs.  In  general,  problems  of  higher  cost  are  more 
difficult  to  solve. 

The  cost  cA  is  given  by  the  equation 

cA  - min  (cfA  + csA  + cdA  , k) 

vhere  A is  the  output  of  an  element  of  type  t and  K is  an  input 

parameter  usually  taken  to  be  32,000.  The  term  cfA  is  that  cost 

contribution  due  to  the  logical  properties  of  t;  csA  is  that  cost 

distribution  due  to  the  side  effects  of  setting  A to  1,  and  cdA  is 

a constant  cost  associated  with  the  element  type  t.  Usually  we  set 
cdA  to  the  values  shown  below: 


A similar  equation  holds  for  cA. 

Calculation  of  cfA  and  cfA 
Gate  functions: 

Consider  a gate  having  input  lines  i ■ l,2,...,n.  Let  the  output 
be  A.  Then  we  have  the  following  results. 
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I 


O 


f 


ABB  gate 


A ■ 1*2  •••  n 


■ i 


A ■ 1 + 5 + ...  + n 


cfA  “ min  {cl} 
i 


BAND  gate 


Interchange  A and  A 


OR  gate 


A **  1+2  + •••+n 


cfA  - min  {ci} 
1 


A - I*  2 


...  n 


NOR  gate 


Interchange  A and  A 


We  represent  a sequence  of  input  vectors  hy  x(l),x(2  ),x(3), 
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Then  the  cost  of  this  sequence  is  the  sum  of  the  cost  of  each  vector. 

Also,  since  C (t ) ■ (C  - 0),(C  » l),  then 

P P P 

, cC  (t)  « cC  + c C 

• P P P 


Flip-flop  (JK  positive  edge  triggered;  S,R  asynchronous): 


cfQ  ■ min  (cS  + cR, 


direct  set 


cJ  + cK  + 2cS  + 2cR  + cC  + cC  , set  F/F 

P P 

cJ  + cQ  + 2cS  + 2cR  + cC  + cC  ) set  or  trigger 

P P 


cfQ  ■ min  (cR  + cS, 


cK  + cJ  + 2cS  + 2cR  + cC  + cC  , 

P V* 

cK  + cQ  + 2cS  + 2cR  + cC  + cC  ) 

P P 


Counter : 


•21 


} 

r 


Q.  - (RL,  RLE  C (t  ) ) 
A p 

+ Rl  P. 


reset  and  increment 


load  a 1 


+ RLE  Cp(t )QA 


Qa  *»  0 and  increment 


reset 


+ RL  P. 


load  a 0 


+ RLE  C (t )QA 
P A 


QA  ■ 1 and  increment 


Therefore 


'I  » 


min  (c£  + 5cL  + 2(cR  + cE)  + cC  + cC  , 

P P 

cR  + cL  + cP.  , 

A 

2(cR  + cL  + cE)  + cC  + cC  + cQ. ) ; 

P P A J 


8 O 


load  a 1 


+ RLE  C (t  )QA 


Q.  * 1 and  increment  ; 


reset 


+ Rl  P„ 


load  a 0 


+ RLE  C (t ) (Q aQb)  Qa  ■ Qg  » 1 and  increment; 


ECPQaWd 

cE  + cCp  + cQa  + cQg  + cC^  + cQp  . 
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Consider  a counter  having  the  input  costs  as  shown  below.  For 
simplicity  we  let  cQ  ■ cfft  for  Q * A and  J. 


(0,«) 


cQa  - min  (l  + 3 • 0 + 2(l  + l)  + 1 + 1 , 

1 + * •••  , 

2(1  + 0 + l)  + 1 + 1 + c§A  ) 

cQa  - min  (7,6+  cQA) 

cQa  ■ min  (l  + 0 , 

1+  »+  ...  t 

2(1+  1+  l)+  1+  1+  cQa) 

■ min  (l , 8 + cQA ) 

■ 1 

cQg  ■ min  (l+  • + •••  , 

2(1 + 0+1) +1+1+  cQa) 

■ 6 + cQa 
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c$g  - min  (l  + 0 , 

1+  «°+  •••  , 

2(l  + 0 + l)  + 1 + 1 + cQA  + cQg) 

- 1 . 

In  summary,  we  have 


cQa  ■ min  (7,6+  c$A) 

(1) 

c«A  - i 

(2) 

1 r 

o«b  - 6 * c«A 

(3) 

and 

eQg  - 1 . 

OO 

Substituting  (2)  into  (l)  we  get 

i o 

CQa  - 7 

(5) 

and  substituting  (5 ) into  (3)  we  obtain 

• 

cQg  - 13  . 

(6) 

To  verify  this  result,  we  see  that  to  set  = 1 we  must  reset 
the  counter  (RL),  and  count  up  two  times,  i.e., 

(RLE  C (t  )f  . 

P 

This  gives  us  a cost  of 


c(RLE  Cp(t)f  - 2(2(1  + 0 + 1)  + 2) 

- 12  , 

hence  a total  cost  of  13. 


Shift  register: 


Legend: 

Device  resets  on  R ■ 0. 

Device  loads  on  L = 0. 

Device  shifts  right  on  S - 1 and  C (t ). 

P 

Data  in  is  D. 

^ - RLS  Cp(t  )Q1_1  shift  in  a 1 

RL  P 


load  a 1 


reset 


> 


i f 


f 


I r' 

l v- 


\ - Rl 

MS  C (t  shift  in  a 0 

Rt  load  a 0 

where 

9o  - D. 

cfSt  - min  (2  (cR  + cL  + cS ) + cCp  + cC^  + cQ^  1 , 

cR  ♦ cL  + ) 

«*I  ■ min  (cR  + cL , 

2 (cR  + cL  + cS ) + cC  + c£  + eft.  , . 

P P i-1 

cR  + cl  + cPt ) . 

Computation  of  Side  Effects  Factor  csA 

When  a line  is  set  to  a 0 or  1,  due  to  its  fan-out  this  line 
setting  may  affect  many  parts  of  a circuit.  We  refer  to  this  phenomenon 
as  side-effects. 

Example: 

Consider  the  circuit  shown  below. 
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r 


o 


To  Justify  F ■ 1,  if  we  set  B ■ 1 we  force  E ■ 1.  But,  if  we 
Justify  F ■ 1 by  setting  C ■ 1,  we  do  not  force  G to  a logic  value. 
Hence,  we  desire  cB  > c C.  That  is,  we  penalize  B because  of  its 
side-effects. 

The  value  of  csA  can  be  computed  as  follows.  Let  line  A have 
logic  value  8 « {0,1} . Simulate  the  circuit  with  the  initial  condi- 
tion A ■ 8 and  all  other  lines  at  x.  The  contributions  to  csA  are 
as  follows: 

1)  a line  B (output  of  a gate)  set  to  0 or  1 contributes  1; 

2 ) a line  B at  x having  m binary  input  has  a contribution 
of  m/n,  where  n is  the  number  of  inputs  to  gate  B; 

% 5)  if  a flip-flop  B is  set,  or  reset,  its  contribution  is  4. 
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gwibcr  EjGi)  of  flip-flops  vhlch  change  state  Is  given  by  the 

9 expression : 


EiM  - 

o z 


o 
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c 


o 


▼ 


Note  that: 

11m  E(n)  - 2 ; 

n -*• 

\ 

R ■ 0:  Resets  all  Flip-Flops 

We  assume  that  half  the  flip-flops  change  state  when  the  counter 
is  reset,  hence  the  cost  contribution  to  csR  is 

!(*♦«> 

e 

where  G Is  the  average  side  effect  associated  with  an  output  line  of 
the  counter,  and  U is  the  cost  associated  with  changing  the  state  of 
a flip-sflop. 

L ■ 0:  fturallel  Load 

We  assume  that  on  the  average,  half  of  the  flip-flops  change  state 
when  a parallel  load  occurs.  The  side-effect's  contribution  to  cL  is 
thus 

-29- 


O 


csL 


- f + G)  . 

E»C_(t)RL:  Increments  a Counter 

p 

When  a counter  Is  Incremented,  E(n)  flip-flops  may  change  state. 
Hence,  the  side-effect's  cost  Is 

* E(n)(4  + G)  . 

Side-Effects  in  an  n-Bit  Shift  Register 

THEOREM  2 : When  an  n-blt  shift  register  is  shifted,  the  expected 

number  Eg(n)  of  flip-flops  which  change  state  is  given  by  the 

expression 

I 

Vn)  " f * 

I The  results  of  'reset'  and  'load'  are  the  same  as  for  the  counter. 

For  the  case  of  a right-shift,  the  side  effect's  cost  is 

, | 0)  . 

Using  these  results,  the  cost  given  earlier  for  the  counters  and 

shift  registers  can  new  be  modified  to  include  the  side  effect's  cost. 

I 

The  side  effects  of  R and  L are  already  in  cR  and  cL,  but 
we  must  add  the  side  effects  due  to  the  increment  process,  since  this 
is  a function  of  both  E and  C (t ).  Hence,  whenever  we  increment  we 
add  the  cost,  2{h  + g),  to  the  cost  given. 
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Hote  that  primary  inputs  only  have  a side  effect's  cost*  Their 
function  cost,  cf , is  0. 

The  coats  for  a circuit  are  computed  twice,  as  explained  below. 

Case  1:  Neglect  side  effects  in  registers  due  to  master  set, 
reset  and  load. 

Case  2:  Use  all  side  effects. 

The  reason  for  employing  these  two  cases  is  the  following.  When  initial- 
ising a circuit  we  do  not  really  care  what  state  the  flip-flops  are  in  as 
long  as  they  are  not  u,  hence  using  master  reset  and  set  is  a good 
strategy,  and  we  should  not  penalize  their  use  because  of  their  large 
side  effect's  cost. 

However,  in  the  middle  of  test  generation,  if  we  need  some  flip- 
flop  to  be  in  the  state  0,  it  could  be  disastrous  to  achieve  this  state 
via  global  reset  line  , since  it  would  also  effect  all  other  flip-flop 
settings.  Hence,  for  this  case,  side  effects  are  used. 

A flow  chart  for  assigning  cA  and  cA  to  all  lines  in  a circuit 
is  given  in  Figure  2 .1,  where  we  assume  the  rules  for  calculating  the  cost3 

for  each  element  type  are  known. 


iitatlon  of  D-Drive  Costs 


The  cost  dA  associated  with  a line  A is  an  estimate  of  the 
difficulty  of  driving  an  error  signal  (d)  on  line  A to  a primary 
output  (po).  The  basic  strategy  for  computing  dA  is  shown  next. 


O 
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FIGURE  2,1  : Computation  of  Line  Costs  cA  and  cX 


O 


1 


We  have  that 

dOUT  *■  min  (dX-^dXg »dXj)  . (l) 

Then 

dA  - dpA  + dOUT  + dt  . (2 ) 

Here  ve  see  that  A is  an  input  to  element  E whose  output  fans -out 
to  three  elements. 

Hie  D drive  costs  associated  with  the  output  line  X from  E 
to  these  three  elements  is  dX^  dXg,  and  dXy  respectively.  Using 
formula  (2  ) above,  then  the  cost  to  drive  a D on  line  A to  an  out- 
put consists  of  three  components,  namely: 

1)  dpA  - the  cost  of  propagating  a D frcro  A to  X,  the 

output  of  Ej 

2 ) dOUT  - the  minimum  cost  of  propagating  the  D at  X to  a po, 
and, 

3)  dt  - a cost  associated  with  the  element  type  t of  E. 
Usually,  dt  increases  with  the  number  of  clock  times 


v 


J 
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required  to  drive  the  D at  A to  X 


1*  Primary  outputs: 


dA  ■ 0 if  A is  a primary  output 


2 • Gate  functions : 


AND  gate  - assume  inputs  1,2 , ,n,  and  a D on  line  i. 


Then 


di  - ci  + dOUT  + 1 . 

J-l , tfr 

This  follows  since  all  other  inputs  to  the  gate  must  be  a 1. 
RAND  gate  - same  result 


OR  gate 


di 


^ cl  + dOUT  + 1 


3.  JK  flip-flop  with  direct  set  and  reset  (negative  logic)  and  positive 
edge  triggering: 

Case  A)  Propagation  of  J ■ D to  output. 

Below  we  see 


-3^ 


i 


I I 


0 


-d 


i 

\f ) 


1 


that  to  propagate  a D (or  5)  from  J to  Q or  5 we  need 
1*  Q ■ 0,  and  a clock  pulse.  Then  if  J « 0 we  have  Q+  =>  0, 
if  J ■ 1 we  have  Q+  **  1.  Therefore, 


dJ  = min  (32000  , cQ  + 2cS  + 2cR  + cC  (t  ) + dOUT  + 2 ) 
Similarly,  we  have 

dK  - min  (32000  , cQ  + 2cS  + 2cR  + cC  + dOUT  +2) 

P 

dS  - min  (32000  , dOUT  + 2 + cR  + cS 

+ min  (cQ  + cC  , 

P 

cQ  + 25p, 
cR  , 

OCp(. ) 

+ min  (c  J + cK  , 
c J + cQ , 
cR  + cK  + cQ))) 


S =*  R » 
and 
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A D on  any  input  line  can  be  propagated  to  any  output 


Hence,  we  can  form  a matrix  of  propagation  costs  shown  below.  The 


entry  shown  is  the  cost  to  propagate  a D on  line  P to  the  out' 


Once  all  entries  in  this  matrix  have  been  calculated,  a 


OUTPUTS 


The  following 
dp(PA 


dp(FA  - dg). 


are  illustrations  of  the  derivation  of  a few  of  these 


To  propagate  a D from  PA  to  QA  we  need  the 
algorithm  LOAD,  l.e. 

C (»)E  R I . 

P 

Hence 

dp(P.  -*  Q.  ) ■>  cC  + c6  + cR  + Ct  . 

A A p 

To  propagate  a D from  PA  to 
algorithm  LOAD  , INCREMENT  , l.e. 


Qg  we  need  the 


C (♦  ) R(EL,  EL) 


Hence, 


dp(PA  -»  %)  - 2cC  (t ) + 2R  + cE  + cL 


+ cE  + cL  • 


dp(L  “*  Qj):  To  propagate  a D from  L to  i e {A,B,C,D}, 

we  must  set  - 8,  P±  - 8,  L » D,  and  clock. 


Cp(t)  E R^  + Pt) 


Hence, 


’ 0 

I 


dp(L  -»  Q±) 


- cCp(t ) + cE  + cR 
+ min  (cQ^  + cl^  , cQi  + cP^ 


dp(E  “*  Qa):  Algorithm  - INCREMENT,  i.e., 


C (t ) "E  - D"  RL  . 


Hence, 


dp(E  -»  Qa)  - oCp(t ) + cR  + cL  . 


Note  that  once  we  have  QA  m D,  we  can  propagate  this  D to 
via  an  INCREMENT  command,  hence: 


dp(QA  - Qg)  - cC  + cE  + cR  + cL 


Similarly, 


dP<S  "*  QC)  ■ cCp  + cE  + cR  + cL  + cQ^ 

*(«C  - V * 00^  + cE  + cR  + cL  + cQ^  + cQg 

dp(Qo  - V ■ cCp  + cE  + cR  + cL  + cQa  + cQg  + cQ^,  • 

These  results  thus  allow  us  to  calculate  the  cost  of 

a)  propagating  a D from  an  input  line  into  the  counter  as  well 

as 

b)  propagating  a D from  one  flip-flop  to  another. 

The  component  dt  of  the  cost  dA  is  given  below: 

S S;  «p  c, 

2 3 4 5 5 

PB  | 5 2 3 4 5 

4 5 2 3 5 

3 4 5 2 5 

22225 
2 3 4 5 5 

2 3 4 5 5 

2 2 2 2 2 

dt 
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n 


5»  Shift  register: 


For  a shift  register  we  propagate  D’s  into  the  register  via  E, 
Cp'  or  R*  We  can  propagate  a D from  to  qi+1  via  a shift 
operation. 

dp(R  -»  Qa)j  Algorithm::  First  set  Q1  - 1 and  then  reset, 

i.e.,  Qj,  - 1 j R . 


dp(R  -»  Q±)  » cQa  + c8  . 


dp(E  -*  0^:  Algorithm::  For  E - D,  clock,  i.e.,  Cp(t ) . R 


dp(E  -*  o ) - 2cR  + cC  + cC 
* P P 


<*p(cp  -»  E):  Algorithm::  Fora  D or  Cp,  enable  shift,  i.e., 
Eli  , hence 


dp(Cp  -*  E)  - 2 (cE  + cfc)  . 


Hate  that  we  can  propagate  a D from  to  Qi+1  via  the 


( 


11  <- 


algorithm, 


E Cp(t ) R , 


hence. 


4p«*i  - Qi+1)  - cCp  + cCp  + 2 (cE  + cR)  . 

When  driving  a D from  E or  Cp  to  we  use  the  dt  cost 
given  by  the  equation  dt  = i.  Driving  a D from  R to  has  a 
dt  value  of  2. 

In  Figure  2.2  we  show  the  flow  chart  for  computing  dA  for  an 
entire  circuit,  where  we  assume  the  rules  for  calculating  the  dA  costs 
for  each  element  are  known. 


O 


Set  all  D drive  costs  to  32000 


I r 


C 


i 


t. 


[SURE  2 .2 : Computation  of  D- Drive  Costs 
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3.  TIMING 


3.1  Basic  Definitions  on  Time 

An  essential  component  in  TEST/80  is  the  way  time  is  handled.  There 
are  three  major  units  of  time  used  in  TEST/80,  namely  a time  frame  (TF), 
a phase  (cp),  and  a period  (p).  In  addition,  the  concept  of  a computation 
step  (s)  is  closely  related  to  time. 

Whenever  a new  input  vector  is  applied  to  the  circuit  a new  time 
frame  (TF)  is  defined.  We  denote  consecutive  time  frames  by  t,  t+  1, 
t+2,...  . Time  frames  which  need  not  necessarily  be  consecutive  are 
denoted  by  t^,tj. 

Example; 

C * t is  equivalent  to 


2 time  frames 


A time  frame  consists  of  one  or  more  phases . For  practical  reasons 

we  limit  the  number  of  phases  in  a single  time  frame  to  1,  2,  or  3* 

The  number  of  phases  in  a TF  is  determined  dynamically,  and  changes  from 
time  frame  to  time  frame.  A line  in  a circuit  can  have  only  one  value 
(static)  during  any  one  phase.  We  denote  the  phases  1,  2,  3 by  <P^,  <P2 , 
Vy  and  <P4(t)  is  phase  1 of  TF  t. 

A phase  is  the  smallest  unit  of  time.  During  a single  phase, 
elements  are  assumed  to  have  O-delay. 
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A 


int,  denoted  by 


> 


is  a device  such  that  if 


\ 


<r 


a signal  enters  it  during  phase  the  signal  exits  at  phase  <Pi+1* 

Every  circuit  must  be  modeled  such  that  each  asynchronous  feedback 
loop  contains  a phase  element.  There  are  two  types  of  phase  elements, 
denoted  by  and  <Pg.  A ^ element  introduces  delay  only.  A <Pg 
element  is  in  a feedback  line  and  represents  a state  variable,  such  as  in 
a latch. 


Example  3.1: 


A B Phase 


1 0 <fL 

0 1 <p3 

1 0 94 


This  situation  (oscillation)  is  not  allowed  for  in  TEST/80  STIM 

GEN. 

The  TEST/80  STIM  GEN  algorithm  allows  for  a maximum  of  three 
phases  per  TF,  hence,  for  two  changes  in  the  value  of  a line  in  response 
to  a primary  input  change.  We  Illustrate  this  "worse"  case  situation  by 
the  circuit  shown  in  Figure  3.1*  Let  the  initial  state  of  the  counter  be 
Ql  ^ Qj  Q4(lsb)  - Hi0  , Cr  - 0,  E-l  and  Cp  - 0.  Let  Cp  now 
go  to  a 1,  hence,  defining  a new  time  frame. 
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(c) 


O 


FIGURE  3.1:  Circuit  Operation  Illustrating  the  Need  for  Three 

Biases. 
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The  following  events  occur: 


9 


t 


i 


ft 


© 


The  clock  causes  the  counter  to  increment  to  state  1111. 

This  forces  * 1,  hence  R ■ 0.  This,  in  turn,  resets 

the  counter  to  the  state  0 0 00  and  hence,  C * 0 and  R 

r 

becomes  a 1. 

Hote  the  need  for  the  phase  element  in  the  feedback  loop. 

In  Figure  3.1(b)  we  show  how  the  signals  may  propagate  through  the 
circuit  as  a function  of  time.  Whenever  a signal  propagates  through  a 
phase  element,  we  have  what  is  known  as  a phase  boundary,  denoted  by  <PB  . 
The  final  steady  state  value  defines  the  phase  boundary  <PF,  which,  in 
turn,  denotes  the  end  of  the  time  frame.  Assume  a phase  element  has  a 
delay  of  A units.  Recall  that  all  other  elements  have  0-delay/phase, 
hence  the  "real  time"  timing  diagram  of  Figure  3.1(b)  can  be  mapped  into 
our  model  "phase"  timing  diagram  of  Figure  3.1(c).  Note  here  that  transi- 
tions during  the  same  phase  occur  in  0 time.  We  see  that  during  the  time 
frame  in  question  changes  value  once,  while  and  R change  twice. 

For  a synchronous  circuit,  there  is  only  one  phase  per  TF,  and 
hence,  we  can  Ignore  the  concept  of  phases  and  deal  only  with  time  frames. 

To  help  clarify  the  modeling  process,  in  Figure  3.2  we  illustrate 
a few  examples  where  no  phase  element  is  needed,  and  in  Figure  3*3  we 
illustrate  a few  examples  where  phase  elements  are  needed.  In  these 
examples,  C represents  an  arbitrary  combinational  circuit. 

Phase  elements  are  also  used  to  introduce  delays,  and  hence,  change 
a critical  race  situation  to  a noncritlcal  race  (see  Figure  3.M* 
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(b) 


A 


(Synchronous 

feedback) 


FIGURE  3.2 : Circuit  Configurations  Where  no  Riase 

Elements  are  Required. 


I 
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A period  la  a unit  of  time  which  corresponds  to  a phase.  Sequences 
in  our  solution  library  are  stored  by  period,  i.e.,  if  X = x^Xg,...  is 
a sequence,  then  x^  corresponds  to  period  p^,  x,  to  period  Pi+1»  x^ 
to  period  etc.  During  test  generation  we  map  periods  into  phases, 

which,  in  turn,  get  mapped  to  time  frames. 

3.2  Problem/Solution-Space 

During  a phase  a line  can  have  only  one  value,  namely,  0,l,u,x,D 
or  D.  Initially,  most  lines  are  assigned  the  value  x at  the  beginning 
of  a phase  and  as  computations  are  carried  out  the  x's  are  assigned 
values  from  the  set  {0,l,u,D,D}. 

In  deriving  a test  of  a fault  in  a circuit,  four  major  operations 
are  carried  out,  namely,  selecting  an  initial  error  generating  stimuli, 
D-drive,  line  Justification  and  implication.  Implication  is  actually 
carried  out  simultaneously  with  the  other  three  concepts. 

For  every  primitive  element  we  have  functional  algorithms  or  tables 
for  carrying  out  these  four  basic  operations.  Except  for  establishing  the 
initial  D in  the  circuit,  test  generation  consists  of  a sequence  of 
solving  two  basic  problems,  namely, 

l)  drive  a 0 through  an  element,  and 

2 ) Justify  a 0 or  1 at  the  output  of  an  element. 

We  denote  the  problem  of  Justifying  line  i at  value  v € {0,1} 
by  the  notation  J : i - v. 

We  denote  the  problem  of  driving  a D (or  D)  from  line  i to  J 
by  the  notation  D : i -*  J. 

If  line  i is  set  to  the  value  v e (0 , 1 , u , D , B}  due  to 
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Implication,  we  denote  this  fact  by  I : 1 *-  v 


We  next  illustrate  some  problem/solution  pairs. 


A.  Example  of  one-period  solution: 


A 

B 


Problem: 


Solution: 


J : C « 0 


0 p. 


B.  Example  of  two- period  solution: 


Problem: 


J : Q 


Solution:  J K 


X X 

1 o 


1 

0 


■i-1 


The  quantities  p^  and  can  be  either  part  of  one  time 

frame,  or  p^_^  can  the  ^ast  P*168®  °f  tinie  frame  t^_^, 
and  p^  the  first  phase  of  tj. 


C.  Example  of  multiple-period  solution: 
COUNTER 
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Problem  : J : ^ ■ 1 , Qg  -1  (assume  initial  state  is  « o) 

Solution:  (E(lx)  C (t  ) , H°)5  (functional  notation) 

Solution  with  no  "holds": 

(E(lx)  Cp(t))3 

Solution  by  periods  (no  holds): 


*1-3 

pi-4 


six  period  solution 


Finally,  it  is  important  to  note  that  a solution  to  a problem  usual- 
ly defines  a new  set  of  problems. 


Example : 


Problem  1 : J : D ■ 1 

Solution  : C ■ 0 


Problem  2 : J : C ■ 0 


Solution  : B ■ 1 


Problem  3:  B * 1 

Solution  : A ■ 0 (solved  since  A is  a primary  output. 


© 


By  transitivity,  solving  Problem  3 in  turn  solves  Problem  2,  which 
in  turn  solves  Problem  1.  (For  simplicity  we  have  ignored  implication.) 

Let  P be  a problem,  e.g.,  a J or  a D.  A computation  steps  s 
consists  of  retrieving  from  the  functional  modeling  library  a prestored 
solution  for  this  problem,  and  applying  this  solution  to  the  circuit. 

This  application  either  leads  to  an  inconsistency,  in  which  case,  back- 
tracking occurs,  or  it  is  consistent.  If  consistent  then  some  new  J 
problems  may  be  created;  new  D- drive  frontiers  exist;  and  some  lines  are 
set  to  0,l,u,D,D  via  implication.  The  total  set  of  new  line  values 
so  defined  are  all  related  to  what  is  known  as  the  step  s computation. 

We  number  our  computation  steps  s ■ 0,1,2,... 

Example: 
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Step 

1 

2 

Line  Number 

3 4 5 

6 

7 

8 

Comment 

1 

X 

X 

D 

X 

X 

X 

X 

X 

Initial  state 

2 

X 

X 

D 

1 

X 

D 

X 

D 

D- Drive  & implication 

3 

X 

X 

D 

1 

1 

D 

D 

D 

D- Drive 

4 

• 

• 

• 

1 

X 

D 

1 

1 

D 

D 

D 

Justify 

I 

1 


1 r 


• Step  2 : Problem  solved  D : 3 -*  6 

Step  3:  Problem  solved  D : 3 -*  7 

Step  4:  Problem  solved  J : 5 * 1 

3*3  Dynamic  Assignment  of  Biases  and  Periods 

When  carrying  out  D- drive  and/or  forward  Implication,  signals  are 
propagated  forward  in  time,  hence,  there  is  no  problem  in  assigning  problem 
solutions  to  phases,  namely  cp„  follows  cp  and  <P,  follows  cp„ . 

^ 1 j 2 

However,  in  line  justification  a sequence  of  line  assignments  for 
periods  pi_^,  pA  may  be  mapped  as  follows : 


O 


© 


© 


*k  - "i  - pi 


Vi  - - Vi 


^ - <P  - p (J -2.3) 
or 

Vl  "Pi-1 


The  test  generation  algorithm  carries  out  this  assignment  process 
automatically.  We  Illustrate  the  need  for  this  type  of  operation  in  the 
next  example. 

Example ; Consider  the  following  circuit  configuration. 
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Row  we  consider  the  following  circuit  configuration,  and  assume 
the  Initial  state  is  Qg  ■ 10  and  we  desire  to  justify  the  state 


The  solution  computation  is  shown  in  Table  3.1.  Ignore  the  time 

and  phase  demarcations. 

To  justify  Qg  Q ^ ■ 01  we  have  the  two  period  solution  shown  in 
lines  1 and  2 , and  the  new  problem  Qg  Q^  ■ 0 0 . To  solve  this  problem  we 
have  the  one  period  solution  R « 0 (line  4).  By  implication,  we  get 
line  3 (forward  in  time)  and  line  3 (backward  in  time).  The  result  from 
line  3 gives  us  the  new  problem  Qg  Q1  ■ 11.  We  solve  this  on  lines  6 
and  7*  We  now  scan  this  solution  from  line  7 up  to  line  1,  assigning 
lines  to  phases  and  time  frasMS.  Whenever  a primary  input  changes  we  get 
a new  TF ; whenever  we  go  up  a line  due  to  implication  (forward  or 


rn 


c 


c 


r, 


0 


Ci 


reverse)  ve  get  a new  phase.  Note  that  the  solution  R - 0 (p^)  gets 

assigned  to  q>2  of  time  frame  tg. 

3.4  Data  Structures 

3.4.1  The  Value  Matrices  V(L,S) 

The  value  of  each  line  In  the  circuit  at  each  step  In  the  computa- 
tion is  stored  In  a set  of  matrices  V(L,S),  where  l ■ 1,2,...,L  is 
the  index  of  a line  and  s **  1,2,..., S is  a step  number.  L Is  the 
largest  line  number  and  S is  the  current  step  being  processed.  Hence, 

V is  of  fixed  width  (l»)  but  of  dynamic  length  (s).  There  is  one  V 
matrix  for  each  unique  phase.  The  V's  are  constructed  dynamically  on 
demand,  l.e.,  whenever  a new  phase  is  created.  Figure  3*3  indicates  the 
general  structure  for  locating  the  value  V(i,s)  « (0 , 1 , D , D , x , u}  of 
line  i at  step  s for  time  frame  t^  and  phase  q^. 

3.4.2  The  Problem/Solution  Tree  (P/S-T) 

Every  solution  to  a problem  can  be  expressed  as  a sequence  of  line 
settings  (j  and  D subproblems)  specified  to  occur  in  a certain  order, 
usually  by  phase  and  time  frame. 

The  general  format  for  a solution  can  be  illustrated  (approximately) 
by  the  data  structure  shown  in  Figure  3*6. 

Initially  the  q?  and  t entries  are  null.  They  are  assigned 
values  dynamically  by  the  test  generation  algorithm. 
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Figure  3.6 


Example: 

The  solution  shown  below: 


J and  D 
line  settings 

• •• 


would  be  denoted  by  the  structure  shewn  in  Figure  3.7. 

Recall  that  a step  consists  of  solving  a D- drive  line  justifi- 
cation problem.  The  effect  of  the  resulting  assignments  set  some  lines 
to  0 and  1 via  Implication.  We  add  these  line  settings  to  the  cor- 
responding data  structure  as  additional  line  settings  at  the  appropriate 
period. 
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Note  that  the  implications  cannot  be  stored  as  part  of  the  function 
problem  solution;  they  are  a function  of  the  circuit  and  not  the  element 
being  processed.  • 

Given  a problem,  we  can  go  to  one  of  our  functional  algorithms  to 
seek  a solution  subtree.  This  subtree  is  then  appended  to  the  problem  tree 
at  the  node  representing  the  problem.  Each  node  in  the  appended  subtree  now 
becomes  a new  problem.  We  illustrate  this  in  Figure  3.8.  The  header  of 
each  subtree  (solution  to  a problem)  is  a step  number. 
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Wote:  If  a solution  subtree  has  a problem  of  the  fora  J : Y - v, 

for  time  frame  tj  and  phase  4^,  but  Y already  has  the  value  v 
for  time  frame  tj  and  phase  then  this  problem  is  deleted  from 

the  problem  subtree. 

The  preceding  description  is  only  illustrative  of  the  concept 
Involved  in  storing  problems/solutions.  The  actual  data  structure 
consists  of  basic  records  having  the  form  shown  in  Figure  3.9. 

The  data  structure  for  storing  the  implications  (i)  created 
by  step  1 is  similar  to  that  used  for  J and  D.  However,  since  I 
does  not  lead  to  any  problems,  we  do  not  require  the  fields,  S,  SF, 

T and  NS. 


f 


O 


used  to  get  next  solution 


****** 


Problem:  Justify  Qg  - 1 (t^ ) 

Solution  No.  1 : Algorithm  (functional) 


(E-lx)  (D-lx  , Cp(t  ) ; (Cp(t)) 

ALGORITHM  (Translated) 


* 

TF  t± 

E D C 

P 

i = n 

X X 1 

- n-1 

110 

■ n - 2 

X X 1 

- n-  3 

o 

X 

H 

In  Figure  J.10  ve  llluatrate  the  data  structure  for  the  solution 
to  this  problem. 

Let  J(s)  he  a list  of  all  (unsolved)  line  Justification  problems 
which  exist  at  step  s.  This  list  consists  of  pointers  to  problems  in  the 


Ve  assume  a synchronous  circuit  where  each  period  is  a unique  time 
frame. 
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FIGURE  3.10 : Data  Structure  of  Problem/Solut ion  in  P/S-Tree 


P/S-tree.  The  list  Is  sorted  in  descending  order  of  cost 


When  going  from  J(s)  to  J(s  + l),  new  entries  cam  occur  due  to 
D- drive  or  carrying  out  line  Justification;  entries  disappear  due  to  line 
Justification  and  implication.  If  S is  the  final  step,  then  J(s)  is 
empty. 

D(s)  is  a list  of  D-drive  problems  existing  at  step  s.  This  list 
consists  of  pointers  to  1>- drive  problems  in  the  P/S-tree.  The  list  is 
sorted  in  ascending  order  of  cost. 

Consider  a problem  created  at  step  J.  Assume,  during  step  S 
(S  > J)  we  attempt  to  solve  this  problem.  If  all  solutions  lead  to  an 
Inconsistency,  then  problem  J has  no  solution  and  all  entries  in  the 
P/S-tree  and  V matrices  corresponding  to  steps  J , J + 1 , . . . , S must 
be  erased. 
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4.  TEST  GENERATION  ALGORITHM 
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In  thin  section  we  will  describe  the  basic  structure  of  the  test 
generation  algorithm  used  in  TEST/80.  In  essence,  it  is  an  extension  of 
the  D- algorithm  and  employs  path  sensitization*  Its  unique  features  are: 

1*  the  use  of  high  level  functional  primitives; 

2.  the  use  of  preprocessing,  namely,  rate  and  cost  analysis; 

3.  the  use  of  a powerful  Initialization  routine; 

4*  the  handling  of  asynchronous  circuits  via  phases  and  time 
frames;  and 

5*  various  miscellaneous  techniques  used  to  reduce  CRJ  time. 

4.1  Data  Structure 

The  data  structure  shown  in  Figure  4.1  is  symbolic,  i.e.,  it  shows 
the  data  representation  of  the  circuit  to  be  tested.  The  structure  is 
based  upon  the  concept  of  a descriptor.  There  is  one  descriptor  for 
every  primitive  element  in  the  circuit. 


1 

; 

s 

V * 

A 

4 

I -.  . 
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The  fields  In  this  descriptor  are  described  next 


T - type  of  element;  T Is  an  integer  and  Is  used  to  Identify  the 
functional  subroutine  associated  with  the  element  type. 

E - unique  Integer  name  for  element  being  defined. 

VALUE -a  vector  having  n components  representing  the  I/O  lines  and 
state  of  the  device ; each  line  can  take  on  the  value  0 , 1 , x , 
D , D or  u,  hence,  there  are  3 bits  per  element. 

The  next  1-1  fields  have  the  format  (j  , address ) vhere  the 
address  points  to  the  J-th  fan-in  of  this  device.  We  assume  1-1  Inputs. 
The  fields  1, . ..,n  are  associated  with  the  outputs  of  the  device.  The 
format  for  these  fields  is  (k , F^ , address).  F^  Is  the  degree  of  fan-out 
of  the  k-th  signal.  If  Ffc  ■ 1 the  address  points  to  the  element  where  k 
fans -out;  if  F^  > 1 the  address  points  to  the  fan-out  list. 

Figure  4.2(a)  indicates  a portion  of  a simple  circuit  and  Figure 
4.2(b)  indicates  the  corresponding  data  structure. 

The  fields  in  the  vector  VALUE  are  ordered. as  follows: 

The  J-th  field,  1 < J < 1 - 1,  corresponds  to  the  J-th  input ; the 

next  fields  (ordered)  correspond  to  the  outputs;  and  the  final  fields 
correspond  to  the  state  variables  (if  any  exist). 

4.2  Functional  Element  Routines 

There  are  five  functional  level  routines  related  to  every  primitive 
element  type.  These  routines  carry  out  the  process  of  initial  error  gener- 
ation, implication,  D- drive,  line  Justification  and  initialisation.  The 


etmcmpt  of  Implication,  D- drive,  and  line  Justification  are  discussed 
1®  our  Report  Ho.  1-76,  "Functional  Level  Modeling  of  Complex  Elements 
in  TEST/80,  and  initialization  is  discussed  in  our  Report  No.  2-76, 
"Initialisation  of  Digital  Logie  Circuits." 

The  routines  are  usually  called  via  a generic  global 
f operation,  e.g.,  D-drive  trc n J input  of  f/f  to  Q output.  The 

tables  respond  with  generic  sequences  for  carrying  out  the  desired 

operation.  The  translator  naps  these  generic  solutions  into  binary 
j strings. 

^•2.1  Initial  Error  Generation 

• 10  Benermte  a test  for  line  i s-a-l(o),  we  place  a 

' D(t)  lnto  the  llna  valAie  table  for  tine  tn  phase  and  store  the 

problem,  "Justify  line  i at  the  value  a"  in  the  P/S-tree,  where 

• a • o(l)  for  the  fault  s-a-l(O). 

^•2.2  Indication 

» Recall  that  every  element  type  is  associated  with  an  implication 

tsble  or  routine.  Whenever  a line  (input  or  output)  of  an  element  is 
changed  frm  x to  v « {0 , 1 , u , D , 5}  we  check  to  see  if  any  other 
» linos,  currently  at  x,  are  forced  to  a value  due  to  this  change. 


f 


x -»  0 


n 

i 

i 


f-i 


When  line  6 changes  from  x to  0,  denoted  by  6 : x 
have,  by  implication  (forward  I,  reverse  I ) 


a)  I 

b)  I 

c)  I 

d)  I 

e)  I 

f)  I 

g)  I 


7 t x -*  1 

8 : x -»  0 

5 : x -*  1 

k : x -*  1 

2 : x -»  1 

1 : x -»  1 


O 


0 


-»  0,  W€ 
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The  process  of  making  all  line  assignments  due  to  Implication  is 
carried  out  by  the  Implication  routine  in  conjunction  with  the  implication 
tables  for  each  element  type. 

This  routine  is  recursive.  It  is  entered  with  an  implication  event 
which  consists  of 

a)  a line  value  v € (0,l,D,D,u}  which  was  formally  an  x ; 

b)  the  name  of  the  line; 

c)  an  element  E;  and 

d)  the  period  p. 

Implication  then  hands  this  data  over  to  the  implication  tables 
(forward  or  backward)  associated  with  element  E.  These  tables  then  return 
with  new  implications  if  they  exist. 

The  implication  events  to  be  processed  are  stored  on  a stack.  They 
can  be  processed  in  any  order.  New  implications  are  placed  on  the  top  of 
the  stack.  They  are  removed  one  at  a time  and  processed.  Implication 
terminates  when  the  stack  is  empty. 

• 

k.3  Implication  Routine 

1.  Pop  implication  event  from  implication  stack.  If  empty  exit, 
else  continue. 

2.  Send  event  to  appropriate  implication  functional  routine. 

3*  Process  response.  Identify  *3-1  lines  changed  from  x due 

to  implication  event. 

U.  For  all  new  line  settings  create  appropriate  implication 
events.  If  event  created  by  forward  implication  on  line  1, 
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then  create  a new  Implication  event  for  all  fan-out  elements 
of  1.  If  event  created  by  backward  implication  on  line  i 
through  element  E,  create  new  implication  events  an  the 
signal  source  for  i and  all  fan-outs  of  1 except  for  E. 
Ait  all  these  new  line  events  onto  the  implication  stack. 

Go  to  1. 

U.h  D- Drive  Routine 

D- drive  is  a simple  routine  which  selects  a D- drive  event  and 
passes  the  event  to  an  element  D- drive  functional  routine.  The  function 
returns  with  the  appropriate  line  settings  required  to  solve  the  D-drive 
event  problem.  The  D-drive  event  consists  of 

a)  the  element  being  processed; 

b)  the  source  of  the  current  D; 

c)  the  target  for  the  D;  and 

d)  the  current  period  and  step  number. 

The  source  can  be  either  an  input  to  E or  one  of  its  internal 
state  variables.  The  target  can  be  either  a state  variable  or  an  output 

line. 

Example 


Right  Shift  Register 


© 
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]>  drive  event: 


a) 

b) 

c) 


d) 

Solution  number 

i: 

D- drive  event: 

a) 

b) 

c) 

d) 

Solution  number  j : 


U1 

I (input) 
(output) 

pU 

(e  * (lx)  cp(t))4 
U1 

% 

% 

p4 

(e  =ix)  Cp(t) 


k.5  Justification  Routine 


Justification  is  a simple  routine  which  selects  a line  justifica- 
tion event  and  passes  the  event  to  an  element  line  justification  function. 
The  function  returns  with  the  appropriate  line  settings  required  to  solve 
the  Justification  event  problem.  A line  justification  event  consists  of 

a)  the  element  being  processed; 

b)  the  line  or  state  variable  to  be  justified; 

c)  the  value  to  be  justified;  and 

d)  the  current  period  and  step  number. 

As  illustrated  previously  in  the  section  on  timing,  for  asynchron- 
ous feedback  loops,  line  justification  requires  dynamic  assignment  of 
line  settings  to  phases. 
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4.6  Initialization  Routine 


This  routine  is  similar  to  D- drive  and  line  Justification.  Its 
function  is  to  obtain  from  the  functional  element  routines  an  input 
sequence  for  initializing  an  element  E,  i.e.,  to  take  the  element  out 
of  the  u state. 


4.7  Consistency  Routine 

This  is  a simple  routine  which,  when  given  a new  line  value 
v(i)  € {0,l,Uj  5,  D]  for  line  i and  its  associated  time  frame  and  cp, 
checks  with  the  V(L,S ) matrices  to  see  that  the  line  does  not  already 
have  a value  w(i)  = [0,1, u,D,  D},  where  v(i)  f-  w(i).  If  this  is  the 

0 

case,  the  new  line  setting  is  said  to  be  inconsistent  with  previous  line 
settings  and  the  algorithm  backtracks. 


4.8  Definition  of  D- Frontier 


The  D- frontier  is  the  set  of  all  elements  having  a D at  a source 
and  an  x at  a target.  The  source  is  either  an  input  line  or  state 
variable;  the  target  is  either  a state  variable  or  an  output  line. 


4.9  Initialization  Module 


Most  test  generation  algorithms  generate  a self- initializing  test 
for  a fault  i.e.,  a test  which  works  no  matter  what  the  initial  state 
of  the  circuit.  In  our  system,  we  use  an  initialization  module  which 
derives  a sequence  R for  driving  the  fault- free  circuit  to  a known  state 
I.  We  assume  R drives  to  the  state  S^.  We  then  derive  tests  for  faults 
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Case  1:  S - Sf.  For  this  case  S/S  is  our  initial  state,  or  goal  state. 


Case  2 : S f Sf  but  Sf  is  completely  specified,  i.e.,  has  no 


For  this  case  S/S^  has  D entries  and  test  generation  is  relatively 

simple. 


Case  3:  S f Sf  and  is  incompletely  specified,  i.e.,  has  u 

entries . 

This  situation  occurs  when  the  fault  f inhibits  the  initialization  of 
some  state  variables.  We  detect  this  type  of  fault  by  the  following 
strategy.  Let  Y = (y1»y2> • • .,yn)  be  the  state  variables  in  the  circuit, 
and  Y^  = (a^s^,, ..  .,an)  be  a state,  i.e.  a^  e {0,1}  for  all  j.  Let 
R be  an  Initialization  sequence  that  drives  the  fault-free  circuit  C^ 
to  the  state  Y^  ■ (a^eu,, . . .,an).  Assume,  due  to  a "reset”  fault  f, 
that  the  circuit  Cf  is  driven  to  the  state  Yg  *=  (u,a^,...,an)  by  R, 

i.e.,  fails  to  initialize. 

/ 

Case  3a:  Assume  the  actual  circuit  Cf  initializes  to  the  state 

vimh- 

We  can  detect  this  situation  by  forming  the  composite  initial  state 
S/Sf  ■ (D,  ag  , a^,  ...  , aR).  Using  this  as  the  initial  state,  we  can 
construct  an  input  sequence  X for  driving  a D to  a po  . 
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Case  3b:  Assume  the  actual  circuit  Cf  initializes  to  the  state 

V 

In  this  case  X nay  not  detect  this  "reset"  fault.  After  RX  has  been 
applied,  assume  the  composite  circuit  CQ/Cf  is  in  the  state  Y^/Y^. 

Let  T be  an  input  sequence  which  drives  C^/cf  from  Y^/Y^  to  a state 

\/\»  y4  “ (*i  > *?.  • bj  > •••  > hj  e {0,1}.  Now,  if  we  apply 

R the  resulting  state  of  Cf  will  be  Y^  - (a^ , 8g  , ...  , an),  and  X 
will  detect  the  fault  (Case  3a).  Hence,  a test  for  a reset  fault  f is 

RXTRX  . 

If  in  trying  to  get  to  Y^/Y^  or  Y^/Y^  we  get  to  a state  Yj/Yj, 

Y^  / Yy  we  again  have  a D and  thus  propagate  it  to  a po. 

4.10  Preprocessing 

There  are  three  major  preprocessing  routines  employed  by  TEST/80, 

namely, 

1.  Rate  analysis ; 

2.  Cost  analysis;  and 

3.  Fault  Collapsing. 

Rate  and  cost  analysis  have  been  discussed  previously.  The 
results  are  stored  in  appropriate  matrices  for  access  by  the  stimulus 
generation  modules. 

As  the  third  phase  of  the  preprocessing  of  a circuit,  faults  are 
collapsed  based  upon  the  classical  concepts  of  fault  equivalence  and 
fault  dominance. 
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4.11  Unsolvable  States 


It  is  often  the  case  that  certain  states  in  a circuit  cannot  he 
reached.  For  example,  a 4-bit  counter  may  be  used  in  a mod  10  counter. 
In  which  case  the  states 

* 

10  10 
10  11 
110  0 
110  1 
1110 
1111 

cannot  be  realized*  Often,  however,  the  test  generation  algorithm  may  try 
to  Justify  one  of  these  states.  Such  a problem  is  unsolvable,  but  may 
require  a great  deal  of.  computation  before  the  algorithm  exhausts  all 
possibilities  and  comes  up  with  this  correct  conclusion.  To  save  CHJ  time 
we  create  a table  called  UNSOLVABLE  which  stores  all  such  problems.  The 
general  form  for  an  entry  in  UNSOLVABLE  is  (S^ , S2 ) where  S1  and  S2 
are  partial  states  (some  X entries)  of  CQ,  and  for  which  there  exists 
no  sequence  for  driving  CQ  from  S^  to  S,,.  For  our  example,  we  would 
write 


and 


(xxxx  , 10 lx) 
(xxxx  , 11  xx)  . 


Entries  are  placed  into  UNSOLVABLE  in  two  ways.  Initially,  the  test 
engineer  can  enter  states  which  are  unsolvable . After  that,  whenever  TEST/ 
80  fails  to  Justify  a state  configuration,  it  caui  enter  this  configuration 
into  the  table. 

i 
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The  table  is  used  as  follows.  When  a new  time  frame  is  being  pro- 
cessed, the  current  state  S of  the  system  is  computed,  and  compared  with 
entries  in  UNSOLVABLE • If  S is  "contained"  in  an  entry  in  this  table, 
then  the  state  is  unsolvable.  (We  say  1 10  x is  contained  in  1 1 x x 
since  1 1 x x covers  1 1 0 x . ) 

4.32  Solvable  States 


Once  the  TEST/80  algorithm  has  computed  a sequence  X.  . for  driv- 
ing  Cq  from  S^  to  S^,  it  stores  this  sequence  in  a table  SOLVE  for 
possible  future  use.  If  this  problem  is  ever  encountered  again,  the 


sequence  X, . can  be  retrieved  and  used, 
“ij 

Entries 


Entries  are  made  into  this  table  automatically  by  the  program.  In 
an  interactive  mode,  the  test  engineer  can  also  make  entries  to  this  ^ble. 
The  format  for  this  table  is  shown  next. 


MttfciCrtr  4^ 
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4.13  Backtrack  State  Event 

There  are  usually  numerous  ways  (solutions)  to  justify  a line  1 
at  a value  v € {0,1} , or  to  D-drive  through  an  element.  If  a selected 
solution  leads  to  an  inconsistency,  one  must  backtrack  to  try  another 
alternative.  The  mechanism  for  backtracking  and  selecting  alternatives 
is  called  a backtract  state  event.  There  is  one  such  program  module  for 
each  element  type.  These  modules  are  part  of  the  functional  element 
type  routines  and  when  given  information  on  the  last  solution,  produce 
the  next  solution. 

There  are  numerous  ways  for  carrying  out  this  process,  and  differ- 
ent element  types  may  employ  different  techniques.  We  will  illustrate 
next  a few  examples  of  how  backtrack  state  events  can  be  handled. 

Example: 
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Justify 


12 


0. 


Solutions : 1 


0 x x 
x 0 x 
X x 0 


1.  Stack  backtract  statement: 
Place  the  three  solutions 


0 

x 

X 


.82- 


X X 
0 X 
X 0 


on  a stack  when  the 


problem  first  occurs.  Whenever  we  need  a new  solution  we  pop 
the  stack.  When  the  stack  is  empty  we  have  exhausted  all 
solutions . 

2.  Counting  backtrack  state  event: 

Solution  i consists  of  placing  a 0 on  the  i-th  input  line 
to  the  gate.  We  enter  the  routine  originally  with  i * 1. 

Each  time  we  backtrack  we  increment  i.  When  i - k,  we  have 
exhausted  all  solutions. 

In  sunmary,  we  associated  with  each  element  type  a backtrack  state 
event  routine  such  that  given  the  last  solution,  the  routine  allows  the 
functional  element  routine  to  produce  the  next  solution,  if  one  exists. 

k.lk  Flow  Charts  for  STIM  GEN 

In  the  flow  charts  presented,  R stands  for  return,  S for  success 
and  F for  failure.  When  we  call  a routine  with  an  argument  R,  the 
routine,  upon  finishing  its  execution,  sets  R to  F or  S depending  on 
the  result  obtained. 

Figure  k.3  shows  the  general  structure  of  TEST/80. 


: Initialization  Module 


f 


I | 
1 


I O 


r 


* 


o 
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The  algorithms  for  this  module  are  discussed  in  our  Report  No.  2-76. 
They  employ  the  line  Justification  and  functional  Initialization  routines 
associated  with  each  element  type.  At  the  conclusion  of  this  routine  all 
initlalizable  circuits  will  be  in  a. given  state.  Faults  which  inhibit  ' $ 

initialization  are  called  reset  faults  and  are  processed  as  discussed  f 
previously. 
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1.  Clear  P/S-tree;  clear  V(L,S)  tables;  S *-  0;  set  Initial 
state  to  IS  (result  of  initialization  module);  TF  «-  1; 


V *-  1. 

2.  Set  up  initial  D on  line  with  fault;  create  D- frontier; 
create  associated  line  justification  problem  « 
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10.  Carry  out  a check  on  rate  analysis  (r) 

R ■ F — (inconsistent)  backtrack  to  7 
R ■ S — continue 

11.  Set  S S+ 1 and  go  to  2 

12.  Is  S - S*  ? 

Ho  : S S - 1 and  go  to  U 

Yes:  R *-  F and  return. 


I 


5.  FAULT  SIMULATION 


5*1  Introduction 

The  function  of  the  fault  simulator  is  to  determine  the  set  of 
faults  detected  by  a given  test  sequence.  The  attributes  of  the  simulator 
for  TEST/80  are  listed  below: 

a)  Concurrent  fault  simulation; 

b)  Functional  modeling; 

c)  Three-valued  simulation; 

d)  Arbitrary  pr operational  delay; 

3)  Input  skew  buffering  (optional); 

f)  Event-directed,  table-driven;  and 

g)  Dynamic  (real-time)  inputs. 

The  rationale  for  using  functional  level  modeling  is  the  same  as  for 
its  use  in  the  stimulus  generation  module,  namely,  efficiency  (reduction  in 
CPU  time)  and  the  fact  that  for  some  element  types,  such  as  LSI's,  logic 
level  descriptions  are  not  available.  The  rationale  for  using  concurrent 
fault  simulation  is  that  it  appears  to  be  the  most  efficient  fault  simula- 
tion technique  computable  with  functional  level  modeling.  This  occurs  be- 
cause the  same  routine  used  to  simulate  an  element  in  the  good  circuit  is 
used  to  simulate  the  same  element  in  the  faulty  circuit. 

To  be  computable  with  our  test  generation  algorithm,  we  can  employ 
functional  level  models  in  our  fault  simulator.  Hence,  even  though 
standard  gate  types  and  flip-flops  are  primitives,  so  also  are  most  MSI's. 
Finally,  some  IflI's  are  primitives,  such  as  RAM's  and  ROM's.  We  refer  to 


The  a Insulator  will  employ  three-valued  logic,  namely  0,  1 and 

u.  The  u has  the  following  meaning: 

a)  steady  state  value  - unknown; 
h)  input  transition  - unknown  ; and 
e)  transient  • race  or  hazard. 

All  S(t^)  must  he  capable  of  handling  this  logic.  For  simple 
elements,  such  as  gates  and  flip-flops,  S(t^)  can  be  Implemented  via 
zoom  tables*  For  more  complex  elements  the  techniques  discussed  in  our 
TKBT/BO  report  on  functional  level  modeling  can  be  used*  (See  also, 
section  on  implication*) 


5.2  Simulation 


The  simulation  system  structure  Is  shewn  In  Figure  5*1 


State  or  Input  Event,  Etc. 


Output  or  State  Event 


SIMULATION 

EVALUATION 


(esem) 


FIGURE  5.1:  Structure  of  Simulation  System 


The  simulator  module  (SM)  passes  the  following  data  to  the  ESEM 


1)  name  of  element; 

2)  type  of  element; 

5)  input  or  state  event  (value  and  line); 

4)  current  total  state  for  this  element;  and 

5)  fault  number  or  good  circuit  event. 

Using  this  data  and  the  corresponding  S(t^),  the  ESEM  executes 
8(t.  ) and  passes  the  following  data  back  to  the  SM. 


l)  output  events,  and 
2 ) the  tine  at  which  they  occur* 

The  nain  modules  in  the  SM  are  shewn  in  Figure  5.2. 


ESEM 

Interface 


Scheduler 
for  Events 

(time  flew 
mechanism) 


Concurrent 
Fault  List 
Processor 


FIGURE  3,2:  Main  Modules  in  Simulator  Module . 


5.3  Data  Structure  for  Circuit 


The  data  structure  used  by  the  simulator  is  a subset  of  that  used  by 
the  stimulus  generation  module,  i*e*,  it  is  based  upon  the  concept  of 
descriptors  associated  with  each  element.  Each  element  type  is  associated 
with  a descriptor  type.  It  may  be  the  case,  however,  that  several  of  the 
elements  modeled  functionally  during  stimulus  generation  are  modeled  at 
the  gate/flip-flop  level  during  simulation.  Ibis  may  be  done  in  order  to 
Obtain  more  detailed  timing  analysis,  or  else  to  treat  internal  IC  faults. 


The  major  components  in  this  data  structure  are: 


a)  the  fan-out  list  associated  with  each  element  output  (the 


O 


r 


fan- In  list  Is  not  explicitly  needed); 

b)  a pointer  to  the  concurrent  fault  list  (CFL)  for  this 
element;  and 

c)  element  Information  consisting  of  l)  type  of  element, 

11 ) delay  parameters* 

The  actual  value  for  the  lines  associated  with  this  element  for 
both  the  good  and  faulty  circuits  is  stored  in  the  concurrent  fault  list 

5*^  . Simulator  Data  Structure 

There  are  six  basic  line  events  which  can  occur,  namely  O/l , 

1/0,  0/u,  l/u,  u/0  and  u/l*  In  general,  we  will  denote  a line  event 

by  e ■ v1/v2 . Because  our  elements  have  multiple  inputs  as  well  as 
multiple  outputs,  events  at  elements  are  actually  vector  events,  e.g., 
(e^eg,...,^)  would  represent  a vector  input  event  to  a n input 
element.  For  simplicity  of  presentation,  we  will  restrict  our  attention 
to  scalar  events,  l.e.,  Just  one  line  event/element  will  be  processed 
at  one  Instance  of  time. 

The  total  state  of  an  element  will  be  denoted  by  (i , S , ©), 
where  I ■ inputs,  S • internal  state,  and  f)  » outputs.  The  current 
total  state  of  a good  element  is  denoted  by  CGTS,  and  that  of  an 
"error  element"  by  CETS,  where  CGTS  / CETS. 

The  concurrent  fault  list  associated  with  an  element  E has  the 


format  shown  in  Figure  5A 


FIGURE  5.U:  Concurrent  Fault  List  Associated 

with  Element  E 

Is  a list  of  faults,  each  of  which  has  the  current  error  total 
state  CETS  i . Rote  that  CGTS  / CETS 1 for  each  i,  and  CGTS  i / 

CETS J for  i i J. 

The  abstract  format  of  the  events  list  which  resides  in  the 
scheduler,  is  shown  in  Figure  5*5*  Here,  t^  < t^+^.  Each  event  in 
list  t^  represents  a signal  change  v^/vg  on  some  line  i in  either 
the  good  and/or  faulty  circuit*  The  actual  structure  of  an  event  entry 
event  J Is  quite  complex. 


© 


Events  List 
Header 


Event 

1 


Event 

2 


* T 


FIGURE  5.5:  Abstract  Structure  of  the  Events  List 


If  In  the  good  circuit  an  event  eQ  « v^/vg  occurs  at  an  input  A 
of  element  E,  the  same  event  can  be  applied  to  all  the  entries  in  the 
concurrent  fault  list  associated  with  E,  with  the  following  exceptions, 
namely,  faulty  circuits  in  which  the  signal  A is  static  at  the  value: 

a)  v^  - these  correspond  to  the  "new  visible  " faults.  They  form 

•add  candidates"  for  the  CFL  of  E. 

b)  Vg  - these  correspond  to  the  "new  invisible"  faults.  They 

form  "delete  candidates"  for  the  CFL  of  E. 

c)  Tj  - these  faults  have  been  and  will  remain  "visible." 

The  following  notation  is  used: 
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«0  - good  event  - v^/vg 

eQ  - 0 no  activity  in  the  good  circuit 


HV 

« new  visibles 

frii 

HI 

• new  invisibles 

»ji 

V 

■ visible  {f^} 

FE 

■ faulty  circuit  event. 

< 


Here,  i la  the  line  on  which  the  event  occurs ; eQ  Is  the  event 
in  the  good  circuit.  The  FCE  list  accumulates  all  the  events  associated 
with  this  line  hut  are  created  by  faults. 

The  following  functions  access  and/or  modify  the  events  list: 


1 - 

1 

If 

i I 


4 
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1.  NEXT- TIME  (time  , nomore) 

If  we  are  currently  working  on  the  header  for  time  t^,  then 

this  function  returns  the  header  for  t.. , for  the  argument 

1+1 

time.  If  no  such  header  exists,  nomore  ■ 1,  else  nomore  * 

0. 

2.  NEXT-EVENT  (event , nomore ) 

This  function  extracts  the  next  event  for  the  current  time 
t^.  If  no  such  event  exists,  then  nomore  » 1. 

5.  SCHEDULE  (time  , event) 

This  function  Inserts  the  new  event,  event , into  the  list 
for  time  - current  time  + delay. 

5*5  Simulator  Operation 

The  basic  structure  of  the  simulator  module  is  shown  in  Figure 
5«6«  The  main  concept  here  is  to  take  events  off  the  events  list,  process 
them,  update  the  CFL,  and  schedule  new  events. 

Notation 


00T8  - current  good  total  state 
NOTE  - next  good  total  state 


▼ 


! 


START 


FIGURE  5,6:  Basle  Simulation  Flowchart. 
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GETS  - current  erroneous  total  state 
NETS  - new  erroneous  total  state 
CES  - current  erroneous  output  value. 


The  program  far  processing  an  event  (IEV)  given  the  current  con- 
current fault  list  (CFL)  in  order  to  product  a new  event  (CEV)  and 
the  new  concurrent  fault  list  (NCFL)  Is  described  next. 

Concurrent  Fault  List  Processor 


i <? 


begin  1 Initialize  CEV  ■ empty,  NCFL  » empty 
if  1 e^  / 0 then  1 

e^  ■ evaluate  (E  , CGS  , eQ ) 
update  CEV*  fQ*  event  ■ e^  "e^  ■ v.j/v^ " 

for  every  f^  e 3EV*V  do 
begin  2 

remove  (f^)  from  CFL 
add  (fk , CFTS ) to  NCFL 

end  2 

for  every  f j e IEV*NI  do 
begin  3 

remove  (f j ) from  CFL 

if  CETS  f NGTS  then  add  (f^  , CETS)  to  NCFL 
if  «q  / 0 and  CEV  ■ v£  then  add  (fj)  to 
CEV.NI 

end  3 


103- 


i 


If 


i 


10 

I 

;i 


\ 


O 


I ■ 

o 
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for  every  e IEV«NV  do 
begin  4 

If  f j « CFL  then 
begin 

remove  (f^)  from  CFL 
add  (f^ , CETS ) to  NCEL 
end 
elae 

add  (t±  , CGTS ) to  NCFL 
if  e^|<  0 ud  QEV  / Vg 

then  add.  (f  ) to  OEV*NV 

end  4 

else  "nothing” 

for  every  t ^ € IEV*DE  do 

begin  5 

if  f^  € CFL  then 
begin  6 

remove  (f^)  from  CFL 
ej  ■ evaluate  (CETS  , e^) 
end  6 

elee  e'  » evaluate  (COTS  , «*) 

If  NETS  f NGTS  then  add  (f  ^ , NETS ) to  NCFL 
if  ej  i 0 then 

if  e^  f e£  then  add  (t * , e^)  to  OEV*DE 
elae  add  (f  ) to  GEV-V 
elae  "ej  - 0"  . 
if  e£  f 0 then 
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o 


o 


o 


Kv' 

ib;: 


if 

'i-'i 

then  add 

<V 

to 

QEV*NV 

if 

t2 

then  add 

<V 

to 

QEV'NI 

► 

s 

pT 

then  add 

<V 

to 

OEV*  V 

end  5 

for  every  X^ 

€ CFL 

do 

begin  6 

e|  - evaluate 

(CETSi 

, eQ)  add 

(xt) 

to 

NCFL 

if 

®i /o 

then 

o 

if  f*  ®o  then  844  (fij » ®£)  to  OEVDE 
else  add  (X1)  to  OEV*V 
else  "ei  - 0" 
if  ®q  / 0 then 


if 

"Ti 

then  add 

<£±) 

to 

QEV»NV 

if 

V2 

“ Vi 

then  add 

U1) 

to 

CEVNI 

if 

v3 

"Vi 

then  add 

to 

OEV*  V 

end  6 


end  1 


5*6  Fault  Insertion 

When  an  output  event  (good  circuit)  occurs  on  a line,  the  associated 
complementary  s-a  - fault  (s  - a - v2 ) entry  is  added  to  the  NCFL,  if  it 
is  not  already  in  the  list,  and  an  entry  made  into  077*DE. 

When  an  input  event  (good  circuit)  is  processed,  then  the  associated 
s-a  - fault  (s  - a - v2 ) becomes  a "new  visible.” 
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When  an  input  faulty  event  (IEV*DE)  occurs  at  E,  we  check  this 


i 


O 


[O 

( 
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event  to  see  if  it  is  compatible  with  E,  e.g.,  the  fault  line  1 
s - a - 8 cannot  place  the  value  B on  line  i. 

In  addition  to  line  s-a  - fault,  we  can  extend  our  faults  to  in- 
clude state  s - a - values  in  functional  models  of  devices  such  as  shift 
registers,  counters,  etc.  " 

5*7  Initialization 

We  assume  the  initial  value  of  all  lines  not  specified  to  be  0 or 
1 by  the  user,  have  the  value  u.  To  initialize  the  circuit  the  fault 
free  simulator  is  run  in  a zero-delay  mode,  i.e.  A 8 0 for  all  elements. 
Only  events  of  the  form  u/0  and  u/l  are  scheduled. 

5.8  Data  for  Primary  Inputs 

We  assume  that  the  data  for  the  primary  inputs  have  been  previously 
generated,  e.g.,  via  the  test  generation  algorithm.  We  also  assume  that 
real-time  inputs  are  possible.  These  inputs  are  assumed  to  be  pre-placed 
into  the  events  scheduler. 

5.9  Input  Skew  Buffering 

If  it  Is  desired  to  process  input  skew,  then  input  events  of  the 
font  v/v  should  be  modified  to  take  the  form  v/u  (A),  u/v  (A). 

9.10  Miscellaneous  Items 

All  classical  aspects  dealing  with  simulation,  such  as  fault  col- 
lapsing, oscillation  control  and  supression,  etc.,  will  be  handled  using 
classical  techniques. 
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